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INTRODUCTION 


The radar altimeters onboard the Seasat and GEOS-3 
satellites were designed for measuring ocean topography and 
monitoring ocean dynamics (Townsend, 1980). During both 
satellite altimeter missions, considerable data were also 
collected over the Greenland and Antarctica Ice sheets. 

The GeoSclence Research Corporation (GSRC) has been 
evaluating the satellite radar altimeter performance over 
the Ice sheets. The purpose of the evaluation Is to provide 
design parameters for future altimeters with improved ice 
sheet measurement capability. To achieve this goal, GSRC 
is utilizing the existing satellite radar altimeter data 
base to determine typical ice sheet surface slopes as well 
as periodicity and amplitude of surface undulations. The 
geometric considerations for satellite* radar altimeter 
measurements over ice sheets are discussed by Brooks, et 
a1. (1978). 

An area of the southwestern Greenland ice sheet was 
selected for detailed topographic analysis; all the available 
Seasat and GEOS-3 altimeter measurements for this area 
have been incorporated into this analysis. 


STUDY AREA DESCRIPTION 


The study area* shown In Figure 1» lies between lati- 
tudes 64.5*N and 65.5”N, and extends from 45.S*W to 48.0*U 
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longitudes. This 13,000 km area of southwestern Green- 
land Is west of the north-south trending Ice sheet crest 
line. The area's mid-latitude approximates the northern- 
most GEOS-3 coverage latitude (65.13*) so that the number 
of GEOS-3 altimeter groundtracks will be maximized. Fourteen 
Seasat and ten GEOS-3 groundtracks traversed the study area. 
Both satellites are In retrograde orbits and overflew the 
area in an east-to-west direction. Uemura's north-to-south 
traverse of Greenland by dogsled. In the Spring and Summer 
of 1978, passed just a few kilometers to the east of the 
study area. He reported percolation facies snow cover for 
the region (Fushimi, et al, 1979). 
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Figure 1. - Location of Study Area 
in Southwestern Greenland. 
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ALTIMETER DATA PR0CESSIW6 
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SURFACE ELEVATIONS 

The data rate from both altimeter sensors was 10 per 
second, providing measurement spacing along each ground- 
track at aoproximately 700 m Intervals. The surface eleva- 
tions (E) referenced to an ellipsoid of a « 6378.137 km 
and f > 1/29S.257 were computed by 

E = H - A ♦ B 

where H Is the satellite height above the ellipsoid from 
orbital computations , A Is the measured altitude, and B 
Is the error In the computation of H. 

The altitudes (A) were recomputed for all the Seasat 
measurements over the Ice sheet. The recomputations were 
accomplished by repositioning the tracking gate at SOX of 
peak power on the surface return waveform ramp Immediately 
preceding the peak power point. Previous experience with 
retracking the Seasat altimeter waveforms over terrain 
features In the United States has resulted In accuracies 
of m when compared with large-scale maps (Brooks, 1981a). 
The requirement for recomputing the Seasat altitudes over 
Ice sheets are discussed by Brooks (1981b). The GEOS-3 
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radar altimeter measurements did not require recomputation 
as Its tracking system was more responsive to the undula- 
tions of the Ice sheet surface. 

The orbital height error (B) Is generally ±5 m. This 
erroi was computed over sea surfaces adjacent to Greenland 
on a pass-by-pass basis by algebraically subtracting the 
altimeter-derived surface elevations from the GEMlOb geold 
heights . 

The subsequent contouring of the Ice surface was done 
hy hand. The contouring Involved no averaging or smoothing; 
every altimeter-derived surface elevation was considered. 


AUTOMATIC GAIN CONTROL (AGC) RENORMALIZATION 


The Seasat AGC values are based on the average power 
in all 60 waveform samples. This calculation of AGC worked 
well over the open ocean, but the waveform mispositloning 
over terrain undulations resulted in erroneous AGC calcula- 
tions. When the surface return waveform ramp was mi sposi tioned 
towards the early gates, the calculated AGC was too high; 
conversely, when the waveform ramp occurred in the later 
gates of the waveform samples, the AGC value was too low. 

For this study, the AGC values were renormalized 
involving the waveform samples in the vicinity of the re- 
positioned tracking gates. Employing an algorithm developed 
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by Htyne (1981, personal comaunlcatlon) , the 60 w«vtfor« 
samples. Si, for each 0.1 second measurement were rescaled 
by g where 

, . 10. <® - “r)/’0- 

0 • AGC value associated with the 
measurement from the onboard 
tracker 

• a reference A&C, assigned a value 
of 35 for this study 

and ■ S^g. 

Next, the 13 rescaled samples In the vicinity of the re- 
positioned tracking gate (I.e., the tracking gate plus the 
6 gates immediately prior to and the 6 gates immediately 
following the tracking gate) were averaged. This average 
value. A, was then used to compute the renormalized AGC 
as 

AGC ■ 10. logio A ♦ Bf 

Where Is a bias, determined to be -»^18.5d6 by comparing 
renormalized AGC values with onboard tracker AGC values over 
open-ocean areas. 
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RESULTS 
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SURFACE ELEVATIONS 

The resultant Ice sheet surface features are depicted 
in Figure 2. The contour Interval In the Figure is 20 m 
with all the elevations referenced to the ellipsoid. 

The 24 altimeter groundtracks are shown by the thin 
lines; the line interruptions indicate areas over which 
the altimeters lost lock. In every instance of loss-of- 
lock, the altimeter waveforms showed evidence of a change 
in the surface ilope. The altimeter did not lose lock at 
every surface undulation, however. These losses-of -lock 
due to surface undulations have also been observed over the 
Antarctic ice sheet and are described by Brooks (1981b). 

The more highly-inclined groundtracks in Figure 2 are 
those of Seasat; the GEOS-3 groundtracks traversed the area 
in a more direct east-to-west pattern, extending northward 
only as far as 65.13®N latitude. The Seasat altimeter main- 
tained lock over the area 82X of the time, while the more 
responsive GEOS-3 altimeter was locked-on 87% of the time 
over the area. Excluding the altimeter drop-outs, a total 
of 3,328 altimeter-derived surface elevations were computed 
for this study area. 

It was observed during the analysis of the altimeter 


7 






measurements that the Ice surface slopes were not continuous; 
flattened areas or terraces appeared quite often along each 
groundtrack. Further, correlations with adjacent or crossing 
groundtracks frequently showed the presence of terraces at 
the same elevations. The terraces were traced from ground* 
track to groundtrack ; their extent (length and width) Is 
shown by the stippled areas In Figure 2 . Some of the terraces, 
e. g., the one at an elevation of 2,424 m, extend across 
the v’ntire study area. An Ice surfac*" profile orthogonal 
to the surface contours between points A and 6 Is shown In 
Figure 3. The average surface slope between points A and B 
Is 4.8 X 10-3. 

The surface elevations of the eastern portion (64. 5* 
65.5®N, 45.5* - 46.5*W) of the study arei* have been plotted 
in three dimensions, courtesy of W. Krabill of NASA; the 
three dimensional plot, shown as Figure 4, more vividly 
displays these terraces. 

The existence of Greenland ice sheet terraces has been 
observed at certain sun angles by scientists on the ice 
sheet (C. Benson, .'981, personal communication), but their 
extent was not previously known. Certainly, it was not 
anticipated that the terraces would persist up to 100 km 
or more in length. It should be noted that the presence of 
thes'' terraces probably would not have been discerned if 
conventional computer-contouring had been utilized; the 
averaging within each grid would have obscured the terraces. 
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Figure 3. - Ice Surface Profile between Points A and B in 
Figure 2. Vertical Exaggeration is 81.4. 







No previous survey elevations from traverses or Geo* 
celver sites are known to exist within the study area. The 
verification of the surface elevations has been performed 
by comparing the surface elevations from Seasat groundtrack 
874 with a Geocelver site prior to the groundtrack entering 
the study area; this groundtrack passed within 750 m of the 
Saddle North Geocelver site, and was the closest of any of 
the groundtracks to a Geocelver site. The geometry and 
surface elevation correlation between the Seasat groundtrack 
and the Geocelver site Is depicted In Figure 5. The Saddle 
North surface elevation Is from Mock (1976), with no estimate 
of accuracy available. The Geocelver site Is described as 
located In the sad(*1e between the northern and southern 
domes of the Ice s.ieet. The maximum positive surface slope 
Is to the west-northwest (Nock, 1976} which Is consistent 
with the elevation comparisons shown In Figure 5. 

RENORMALIZED AGC 

The renormalized AGC values have been correlated with 
the surface elevations resulting from the waveform retracking. 
Generally, wherever the surface slope Increased, the re- 
normalized AGC values decreased; and wherever the surface 
slope decreased at one of the terraces, the AGC increased. 

As an example. Figure 6 contains a detailed portion of the 
surface elevations from Seasat orbit 806. This Figure also 
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Figure 6. - Example of Correlation between Surface Slope 


and Renormalized AGC. 
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Illustrates the corresponding renormalized AGC values and 
the excellent agreement between surface slope and signal 
level . 
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ALTIMETER PERFORHANCE 


The Seasat and GEOS-3 aU1fflett*s performed quite well 
over this Ice sheet area. Presuming that subsequent satellite 
altimeters to be launched in this decade will be of the 
Seasat-type, the analysis of performance has concentrated 
on the Seasat altimeter. 

A typical waveform from the study area surface as 
recorded by Seasat Is shown In Figure 7. The waveform rises 
abruptly and significant power levels are maintained for 
about 30 gates. Indicating a highly reflective nadir target 
as well as numerous off-nadir surface facets oriented In the 
direction of the altimeter's receiving antenna. This wave- 
form Is typical, also, In that the waveform rise Is mls- 
posltloned by the onboard tracker with respect to the tracking 
( or 0th ) gate. 

The average range velocities for the study area, based 
on the surface contours and the groundtrack orientation with 
respect to the contours, varied from 0-40m/sec with Instan- 
taneous velocities frequently exceeding 100m/sec. An example 
of calculated instantaneous range velocities Is illustrated 
In Figure 8 where, for a portion of Seasat orbit 806, the 
average range velocity approximates -^20m/sec. The corres- 
ponding surface elevations are In Figure 6. 

The Instantaneous range accelerations, again based on 
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the surface elevations, frequently exceeded lOOm/sec 
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Figure 7. - Typical Seasat Altimeter Surface Return Waveform 


for Study Area. 
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Calculated Instantaneous Range Velocities for the Surface Elevations In Figure 



at the surface terraces. The range accelerations for the 
same portion of orbit 806 are depicted In Figure 9. 

As discussed by Brooks (1981b)» the altimeter losses- 
of-lock are not the result of high Instantaneous surface 
accelerations, but are produced by a persistent surface 
change which causes the altimeter waveform to "walk-out" 
of the range gates. This "walking-out" occurs during 0.5 - 
1.0 second time spans. A faster tracker update would solve 
this problem. As long as the waveform stays within the 
range gates, the waveform can be retracked to produce accurate 
surface elevations. 
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Figure 9. - Calculated Instantaneous Range Accelerations for the Surface Elevations In 
Figure 6. 



SUMMARY 


Fourtten Seasat and ten GEOS-3 satellite radar altliaeter 
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groundtracks across a 13,000 km Ice sheet area of south- 
western Greenland yielded 3,328 surface elevations. The 
surface elevations derived from Seasat were recalculated 
based on a waveform retracking algorithm. Analysis of the 
elevations reveals the existence of surface terraces, some 
greater than 100 km In length. Renormalized Seasat A6C 
values were shown to be correlated with surface slopes. 

The Ice surface undulations caused frequent altimeter 
losses-of-lock. A faster servo response rate for future 
satellite altimeters would aid Ice sheet mapping In the 
next decade; as long as the waveforms are retained In the 
range gates, the waveforms may be subsequently retracked to 
achieve accurate surface elevations. 
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